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The prospects of extracting transport coefficients from euclidean lattice simulations are discussed. Some general 
comments on the reconstruction of spectral functions using the Maximal Entropy Method are given as well. 



1. In field theory transport coefficients are pro- 
portional to the slope of appropriate spectral 
functions at zero frequency and zero spatial mo- 
mentum (Kubo relation) . Examples are the elec- 
trical conductivity, 
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where p cm {x,y) = ([j l cm (x) , j l eia (y)]) is the spec- 
tral function associated with the electromagnetic 
current j* m = V^*^, an d the shear viscosity, 



(2) 



where p^(x,y) = ([^(x), 7r«(y)]) with tt m the 
traceless part of the spatial energy-momentum 
tensor. 

The euclidean two-point function and the spec- 
tral function (both at zero momentum) are re- 
lated via a dispersion relation and 



G(r) = 
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with the kernel 



K(t,u) 



r n B (u) + e-" T [l + n B (uj)], (4) 



where ns(w) = l/[e"/ T — 1] is the Bose distri- 
bution. The first attempt to compute transport 
coefficients on the lattice using Eq. (0) was made 
some time ago by Karsch and Wyld M , by fitting 
a three-parameter ansatz for the spectral func- 
tion to the euclidean lattice data. This approach 
was pursued more recently in Ref . |j| . A modern 
way to attack this problem would of course be to 
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use the Maximal Entropy Method Qj. Once the 
spectral function is reconstructed for all to, the 
transport coefficient is in principle determined. 

Two obvious questions are: What is the spec- 
tral function expected to look like at high tem- 
perature? How does the transport coefficient, or 
in general the low-frequency part of the spectral 
function, show up in the euclidean correlator? 

2. To answer the first question, we calculated 
the spectral function relevant for the shear viscos- 
ity at high temperature in weakly-coupled scalar 
and nonabelian theories jy. The results for the 
scalar theory are sketched in Fig. [l]. The contribu- 
tion at higher frequencies, i.e. w > 2m where m is 
the thermal mass, arises from (inverse) decay pro- 
cesses. At large enough u> the spectral function 
increases as p wv (ij) = {oj 4 /96tt) [1 + 2ns(w/2)]. 
The origin of the contribution at lower frequen- 
cies is scattering of fast-moving particles with soft 
bosons in the plasma. The approximate shape of 
the spectral function in this region is given by 
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(w«T), 
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where 7 is the thermal damping rate. The par- 
ticular form in Eq. (^|) is due to the presence of 
poles in the complex energy-plane that pinch the 
real axis from above and below. The denomina- 
tor indicates the distance between these so-called 
pinching poles. The viscosity is determined by 
the slope at zero frequency and is proportional 
to I/7. A quantitative analytical calculation in 
this regime is complicated since due to the pinch- 
ing poles the loop expansion breaks down and all 
ladder diagrams with uncrossed rungs contribute 
at leading order in the coupling Q . In the scalar 
theory, the low- and the high-frequency contribu- 
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Figure 1. Sketch of the spectral function p^(ui) /T A versus u> in a weakly-coupled scalar A</> 4 theory (not 
to scale). The shear viscosity is given by the slope at zero frequency. See the main text for discussion. 



tion match parametrically around the threshold 
for decay u> = 2m. A general ansatz describing 
the low-frequency contribution is 
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where x — lu/T and bi — Ci = 0, i > n for given n. 
For gauge theories the overall characteristic shape 
remains as in Fig. [H but the dependence on the 
coupling constant g ~ A differs somewhat Q . 

3. We now address the second question raised 
in the beginning. At leading order in the cou- 
pling constant the euclidean correlator can be 
computed exactly using Eq. (0), and we find 



Gtttt (t) = ^ { (n — u) [11 cos u + cos 3it] 
3 sin u y- ' 



3 sin 5 u 
6 sinu + 2 sin 3u| 
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with u = 1~ktT . For the scalar theory cii = 
a% = 1, while for an SU(N C ) theory a\ = 602 = 
12(A r ( ? — 1). The r-dependent terms arise entirely 
from the high-frequency part of the spectral func- 
tion. The r-independent term originates from 
the low- frequency part. The fact that the low- 
frequency part of the spectral function leads to a 
constant term in the euclidean correlator is easily 
understood. Since for small frequencies lu <C T 



the kernel (Q) can be expanded as 



(8) 



and the dominating term at small to is r- 
independent, the low- frequency part of the spec- 
tral function corresponds in the euclidean cor- 
relator to a constant term proportional to 
J diu p n7r (u))/uj. This particular constant cannot 
be easily disentangled from the high-frequency 
contribution. This is illustrated in Fig. ^, where 
we plot the analytical result for G 7r7r (r) in the 
case of SU(3). The tiny difference between the 
full and the dashed lines is due to the constant 
term in Eq. (0). Once again, this term originates 
from the low-frequency part of the spectral func- 
tion and reflects J du> p^iui) / u> , not /9 W7r (w) itself. 
Therefore, we find that G ffW (r) is remarkably in- 
sensitive to details of p^i^Lo) when uj <C T and we 
conclude that it is extremely difficult to extract 
transport coefficients in weakly-coupled field the- 
ories from the euclidean lattice. 

4. The results described above are generic and 
not specific for the correlators we considered. The 
fact that the low-frequency part of a spectral 
function corresponds to a constant term in the 
euclidean correlator relies solely on the expan- 
sion (||). In our opinion this result is a poten- 
tial problem for the Maximal Entropy Method 
when the reconstruction of low-frequency parts 
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Figure 2. Logarithm of G^^M/T 5 versus tT in 
weakly-coupled SU(3) (Eq. (fjj), full line). The 
dashed line is the result without the constant 
contribution originating from the low-frequency 
part of the spectral function. The inset shows a 
blowup around tT = 0.5. 



of spectral functions is attempted. We empha- 
size that the difficulty is not a numerical one, 
e.g. due to a finite number of lattice points in 
the imaginary-time direction. Furthermore, it is 
easy to see that at finite temperature the pres- 
ence of pinching poles in correlators of composite 
operators that are bilinear in the fields is quite 
common. From general considerations it follows 
that the effect of pinching poles is a low-frequency 
contribution as in Eqs. ^]) and Fig. [l]. 

Consider for instance the electromagnetic 
current-current correlator (or any two-point 
function of fcrmion bilinears) in the decon- 
fincd quark-gluon plasma. For large frequen- 
cies a perturbative calculation gives p em (co) ~ 
lu 2 [1 - 2n F (w/2)], where n F {uo) = l/[e w / T + 1] 
is the Fermi distribution and the proportional- 
ity factor depends on the number of fermions, 
colours, etc. For this reason it is customary 
to present this kind of spectral functions as 
p em (tu) / uj 2 . However, due to pinching poles also 
this spectral function is expected to have a struc- 
ture at small frequencies as in Eqs. (||, ^ and 
Fig. |l|, where in this case 7 is the fermion damp- 
ing rate. For very small frequencies the spectral 



function can be expanded as 

Pem(w) 
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where d\ is proportional to the electrical conduc- 
tivity. The coefficients di are nonzero at finite 
temperature only and represent repeated scatter- 
ing of the fast-moving on-shell fermions with soft 
gauge bosons in the deconfined phase of the high- 
temperature plasma. Because of this behaviour, 
p em (uj)/uj 2 diverges as 1/w when uj —> 0. Singu- 
lar behaviour at very small frequencies of spectral 
functions of fermion bilinears that are normalized 
with 1/w 2 has indeed been found || although the 
statistical significance of these results is still un- 
certain. 
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